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Systems of free fermions are classified by symmetry, space dimensionality, and topological prop¬ 
erties described by K-homology. Those systems belonging to different classes are inequivalent. In 
contrast, we show that by taking a many-body/Fock space viewpoint it becomes possible to estab¬ 
lish equivalences of topological insulators and superconductors in terms of duality transformations. 
These mappings connect topologically inequivalent systems of fermions, jumping across entries in 
existent classification tables, because of the phenomenon of symmetry transmutation by which a 
symmetry and its dual partner have identical algebraic properties but very different physical inter¬ 
pretations. To constrain our study to established classification tables, we define and characterize 
mathematically Gaussian dualities as dualities mapping free fermions to free fermions (and inter¬ 
acting to interacting). By introducing a large, flexible class of Gaussian dualities we show that any 
insulator is dual to a superconductor, and that fermionic edge modes are dual to Majorana edge 
modes, that is, the Gaussian dualities of this paper preserve the bulk-boundary correspondence. 
Transmutation of relevant symmetries, particle number, translation, and time reversal is also inves¬ 
tigated in detail. As illustrative examples, we show the duality equivalence of the dimerized Peierls 
chain and the Majorana chain of Kitaev, and a two-dimensional Kekule-type topological insulator, 
including graphene as a special instance in coupling space, dual to a p-wave superconductor. Since 
Gaussian duality transformations are also valid for interacting systems we briefly discuss some such 
applications. 


I. INTRODUCTION 

In this paper we establish equivalences of topolog¬ 
ically non-trivial insulators^^^ and superconductors. 

By means of duality transformations, we show that any 
insulator has a dual superconducting partner, and the 
partners are either both topologically trivial or non¬ 
trivial. We will focus on non-interacting dual partners, 
since general classification schemes exist for free-fermion 
systems.®^® As it turns out, the duality transformations 
of this paper connect systems that are inequivalent from 
the point of view of these topological classifications. This 
is only possible because of the phenomenon of symmetry 
transmutation, by which a duality transformation maps 
a symmetry of one system with one physical interpreta¬ 
tion, say particle number or time reversal, to a symmetry 
of the dual system with a different interpretation. 

From an electromagnetic response viewpoint, insulat¬ 
ing and superconducting phases of electron systems are 
dramatically different. While the insulating phase is 
characterized by a vanishing current-carrying state at 
zero temperature, the superconducting phase supports 
a supercurrent and displays a perfect diamagnetic re¬ 
sponse, the Meissner effect. Yet there is a basic sense in 
which both states of matter are equivalent, since many 
of their defining properties stem from a common factor, 
that is, the existence of a gap in the bulk energy spec¬ 
trum of fermionic quasiparticles. The additional presence 
of gapless, symmetry-protected, extended surface excita¬ 
tions defines operationally their topologically non-trivial 
character. One of the objectives of topological band the¬ 
ory is to classify, based on a few preferred (discrete) sym¬ 
metries and space dimensionality, topologically distinct 
non-interacting (single-particle) Hamiltonians and their 


concomitant gapless edge excitations. 

For systems without gauge symmetries, duality trans¬ 
formations are implemented by unitary mappings, 
and so they preserve symmetries; the symmetries of a 
system are in one-to-one correspondence with the sym¬ 
metries of its dual partner. However, the physical in¬ 
terpretation of a symmetry and its dual image can be 
markedly different. Holographic symmetries'^ constitute 
a most extreme example. For some pairs of dual part¬ 
ners, one of the systems displays boundary symmetries, 
mapped to global symmetries of its dual partner. In this 
case we call the boundary symmetry holographic. This 
phenomenon is remarkable because a symmetry that is 
formally lost in the thermodynamic limit, the holographic 
symmetry, is mapped by duality onto a symmetry that 
may become spontaneously broken in that limit. Because 
of this, not uncommon, example of symmetry transmu¬ 
tation, it is conceivable that a duality may map a parti¬ 
cle conserving system to a non-conserving one, simply by 
mapping the C/(l) symmetry of particle number to a dual 
U{1) symmetry that does not have that interpretation. 

While these arguments are encouraging in the search 
for equivalences of insulators and superconductors, 
there are at least two other obstacles besides parti¬ 
cle (non)conservation. First, in general, dualities for 
fermions will not preserve the quadratic (or Gaussian) 
character of a model system, often mapping free-fermion 
systems to interacting ones. Second, for topologically 
non-trivial systems, even if one were to find dualities 
matching non-interacting dual partners, there is in gen¬ 
eral no reason to expect that these dualities should 
also preserve the locality properties of the quasiparticle 
modes. For example, at zero energy, modes localized at 
a boundary may also be interpreted as boundary sym- 
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metries. Extrapolating from the experience with holo¬ 
graphic symmetries, one would expect these modes to 
become delocalized after a duality transformation. 

As it turns out, both obstacles may be overcome, and, 
as a consequence, there is no fundamental obstruction 
to the construction of equivalences of topological insula¬ 
tors and superconductors in terms of dualities. Section 

II introduces the special class of duality maps that es¬ 
tablishes those equivalences. The starting point is the 
characterization of duality transformations that preserve 
the quadratic fermionic nature of a given model system; 
we will call this transformation Gaussian duality. Next 
we will construct a large class of such dualities in any 
number of spatial dimensions, in order to create a toolkit 
for generating a topological superconductor from any 
given topological insulator in a systematic fashion. In 
other words, given a topological insulator one can always 
find at least one dual topological superconductor asso¬ 
ciated to it. This process is, of course, reversible since 
the duality transformation is an isometry. Hence our re¬ 
sults strongly suggest that there may exist equivalences 
of topological insulators and superconductors across all 
entries of the topological classification table, at least for 
constant space dimension. Dimensional reduction by du¬ 
alities is possible,but we will not obtain any Gaussian 
instance of this phenomenon in this paper. Section II 
ends with the fundamental concept of symmetry trans¬ 
mutation as applied to fermion parity, translation and 
time-reversal symmetries. In particular, we will find a 
quantitative connection between changes in translation 
symmetry and breaking of particle conservation. 

Particular and emblematic examples include the proof 
that the insulating dimerized Peierls^"^’^® and supercon¬ 
ducting Kitaev (at vanishing chemical potential) 
chains are dual partners, and the equivalence of graphene 
to a popular example^®’^® of a weak^^ topological super¬ 
conductor in two spatial dimensions. It is in Section 

III that we present these two prototypical equivalences. 
These dual partners do not simply resemble each other, 
but are isospectral from a many-body standpoint, for fi¬ 
nite lattices and various boundary conditions. No doubt, 
this fact seems odd at first sight, since the spinfull Peierls 
chain for example partially breaks translation symmetry 
but not (the standard) time reversal or particle conser¬ 
vation, while its dual partner, Kitaev’s Majorana chain, 
breaks time reversal and particle conservation but not 
translation symmetry. (In Appendix A we derive the su¬ 
perconducting dual of the m-merized Peierls chain, and 
discuss the differences between m odd and even.) Simi¬ 
larly, graphene displays the symmetry of the honeycomb 
lattice^*^ while it dual superconducting partner sits on a 
square lattice. Remarkably, the Gaussian duality allows 
us to qualitatively understand the difference between zig¬ 
zag and armchair terminations in graphene.The expla¬ 
nation to all these seemingly paradoxical observations is 
symmetry transmutation. In Appendix B we describe a 
simple Gaussian duality mapping a (s-wave) BGS super¬ 
conductor to an insulator in any number of dimensions. 


Another important issue addressed in Section III is 
the locality character of our Gaussian dualities, i.e., the 
problem of showing that localized zero-energy modes 
are mapped to dual zero-modes that are also local¬ 
ized. It is remarkable to have the possibility to gener¬ 
ate localization-preserving Gaussian dualities. In other 
words, there are no holographic symmetries associated to 
the Gaussian dualities of this paper: global symmetries 
map to global symmetries, and the localization properties 
of energy modes are also preserved, though edge modes 
may be shuffled among boundaries. In particular, the 
zig-zag boundary of graphene is exactly dual to the “Ki¬ 
taev edge” of Refs. 18 and 19. We show how topological 
defects and edge states map, and also how the nature of 
those excitations transmutes from (canonical) fermionic 
to Majorana character by duality. Interestingly, we ana¬ 
lytically construct exact (as opposed to asymptotic) zero- 
energy modes for any finite length Kitaev wire when the 
length is an odd number of lattice constants. 

An interesting outcome of our investigation is further 
confirmation that non-trivial topological quantum order 
is a property of a manifold of states interpreted relative 
to a given language (a set of preferred observables^^); and 
not a property of the energy spectrum alone^^ and some 
Hamiltonian singling out those states as energy eigen¬ 
states. We also investigate in Section HI the interplay 
between dualities and topological invariants of the single¬ 
particle Hamiltonian. Indeed, our equivalences are dual¬ 
ity mappings and hence necessarily isospectral.^^ How¬ 
ever, some of these duality mappings connect systems 
with ground states characterized by different topological 
quantum numbers, thus belonging to different topolog¬ 
ical classes. For instance, Kitaev wire model belongs 
to the class D of the Dyson-Altland-Zirnbauer tenfold¬ 
way classification,® while its dual, the dimerized (spin¬ 
less) Peierls chain, belongs to the class AHI. In the past, 
we have studied dualities mapping systems with topo¬ 
logically quantum-ordered ground-state manifolds to sys¬ 
tems characterized by local (Landau) orders. 

In order to take advantage of duality transformations, 
it is crucial to recognize that the many-body, and not the 
single-particle, representation of the system is the rele¬ 
vant one. There is absolutely no doubt that, from a com¬ 
putational standpoint, the single-particle representation 
(e.g., the Bogoliubov-de Gennes equations) is the appro¬ 
priate methodology to adopt in the non-interacting or 
mean-field case. Computationally, it reduces a problem 
of exponential complexity into one of polynomial com¬ 
plexity, thus allowing diagonalization of quite large sys¬ 
tem sizes. However, care must be exercised at the mo¬ 
ment of analyzing properties like particle conservation 
that involve the whole many-body system. In partic¬ 
ular, topological classification schemes and counting of 
many-body zero-modes relate directly to the many-body 
ground-state manifold. These cautionary remarks are 
entirely appropriate since Gaussian dualities connecting 
topologically non-trivial dual partners often seem at odds 
with one form or another of standard wisdom. However, 
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they are entirely natural if one adopts the many-body 
(Fock-space) language of second quantization, and not 
the vector bundle analysis of single-particle Bogoliubov- 
de Gennes Hamiltonian matrices. 

Technically, Majorana operators, defined (up to nor¬ 
malization) as the real and imaginary parts of the canon¬ 
ical fermionic field, generate a complex Clifford alge¬ 
bra naturally represented in Fock space, and our Gaus¬ 
sian dualities are characterized most naturally as isomor¬ 
phisms of these Clifford algebras. The effect in single¬ 
particle (mode) space, where a different, exponentially 
smaller Clifford algebra emerges®’^ is induced a poste¬ 
riori. Crucially, it follows that our Gaussian dualities 
can also be used for investigating interacting many-body 
systems (see for example Section III A 3).^® Several dif¬ 
ferent mathematical simplifications arise when Gaussian 
dualities are investigated in terms of Majorana opera¬ 
tors. These simplify not only the search for equivalences, 
but also the analysis of symmetries, topological invari¬ 
ants and their transformation, and most importantly the 
mapping of boundary excitations. 

At this point it becomes natural to ask about the ex¬ 
tension of our work to bosons, since it is clear that the 
notion of Gaussian duality applies to canonical bosons 
just as well. However, the real and imaginary parts of 
the bosonic field satisfy the Heisenberg commutation re¬ 
lation, and so the theory of Gaussian dualities for bosons 
is bound to be markedly different from that of fermions. 
Due to this crucial technical difference, we defer the sys¬ 
tematic study of bosonic Gaussian dualities to future re¬ 
search. Nonetheless, we would still like to illustrate ex¬ 
plicitly the point that symmetry transmutation is also 
operative in bosonic systems. Hence, in Appendix G we 
describe a duality mapping of a bosonic Mott insulator 
to a quartet superconductor. We also comment briefly 
on the relevance of this example for cold atoms. 

Section IV concludes with a summary and outlook. 


II. GAUSSIAN DUALITIES 


In statistical mechanics, duality mappings are unitary 
transformations that respect the locality structure of the 
many-body Hamiltonian or transfer matrix.In what 
follows we will set up the foundations to establish equiv¬ 
alences via dualities. Particularly, we will character¬ 
ize operator maps relating Hermitian quadratic forms of 
fermions, i.e., Gaussian dualities, including the connec¬ 
tion between such many-body dualities and the associ¬ 
ated transformation of the single-particle Hamiltonian. 
Next we will introduce general techniques to decompose 
a very large class of free fermion models into sums of com¬ 
muting Hamiltonians, and finally we will use this tech¬ 
nique to construct a general class of Gaussian dualities. 


A. What constitutes a Gaussian duality ? 


We are interested in establishing the conditions under 
which a duality transformation becomes Gaussian. Let 
us focus for simplicity on systems of free fermions de¬ 
fined on a lattice. Then, the most general free fermion 
Hamiltonian is of the form 


H = 


^ 1 1 

-I- -I- -Aj^-c^Cj 


( 1 ) 


with one-body and pairing interaction matrices 

K‘^ =K, A^ = -A, (2) 

where ^ is the adjoint, ^ the transpose, and * complex 
conjugation of a matrix. The creation (annihilation) op¬ 
erator of a fermion cj (c^-) in the single-particle orbital 

4>j, satisfying { is labelled by the generic 

subindex j encoding arbitrary quantum numbers, includ¬ 
ing position, spin, orbital/band, angular momentum, etc. 
The total number of single-particle orbitals is L. 
Equivalently, one can re-write H in Nambu form 

1 at/iBdG a-I-^Tr AT, (3) 

where the column vector of fermion operators is given by 

a = , with aj = Cj , ar+j = cj , j = I, • • • , L, (4) 

and the Bogoliubov-de Gennes single-particle Hamilto¬ 
nian (2L X 2L matrix) 

hB^G = = ( 5 ) 

il G ^{K) + \t^ (g) 9(A) -F ir^ (g) 3?(A) +t^ ® 3?(A), 

where r^, v = x,y,z, are Pauli matrices, and 5ft(-)(9(-)) 
denotes the real (imaginary) part of the matrix. No mat¬ 
ter what the specific matrices K and A are, the single¬ 
particle Hamiltonian ft-BdC always anticommutes with the 
antiunitary (particle-hole) operator 

C = /Cr^(g)l, C^ = l, (6) 

i.e., {ft-BdCjC} = 0, where /C denotes complex conjuga¬ 
tion. That means that the single particle energy spec¬ 
trum is antisymmetric with respect to its zero value, i.e., 
a particle-hole symmetric spectrum. By contrast, a chi¬ 
ral symmetry Z^chirai is a unitary transformation that anti¬ 
commutes with /iBdG- For example, if 9(Ar) = 0 = 9(A), 
then Z^chiral = (g) 1. 

Suppose now that the unitary transformation Z^d im¬ 
plements a duality transformation, 

H^^U.HUl ( 7 ) 

meaning that it transforms a local non-interacting Hamil¬ 
tonian H into another, dual , that preserves the prop¬ 
erty of being also local. The map, however, could gen¬ 
erate fermionic density-density interactions for instance. 
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What are the general conditions under which is also 
an Hermitian quadratic form of fermions? 

To answer this question we will recast Hamiltonian H 
of Eq. (1) in terms of Majorana operators 

72i-l = Cj + c], n 2 j = Cj - c], (8) 

such that 7 ^ = 1, r = 1, • • • , 2L. (The notation 

aj = Cj+c], \bj=Cj-c]. ( 9 ) 

will be favored in later sections). Then 

i 1 

H = hrslrls + 2 "’"''^ (10) 

r,s—l 

becomes a quadratic form of Majorana fermions, with a 
2L X 2L matrix h that is real and antisymmetric. 

We can now investigate the dual Hamiltonian. Re¬ 
member that we want the duality map to be Gaussian, 
i.e., should also be a quadratic form of Majorana 
fermions. A naive first, and trivial, attempt would be to 
keep the localization properties identical, i.e. 

. 2L 

^ + -Tr K (11) 

r,s=l 

where the dual operators are related to the originals as 

( 12 ) 

and 7,P is a Majorana fermion operator. This extreme 
local map, although Gaussian, is very restrictive and will 
not allow us to establish interesting equivalences between 
insulators and superconductors. We would like to relax 
the extreme locality constraint and allow for changes in 
the range of the dual matrix. In other words, we would 
like to realize a more general Gaussian duality 

. 2L 

= f^?slrls + ^TrK, (13) 

r,s—l 

for some new Majorana operator 7^ significantly different 
from the dual Majorana 7,? and the original one 7^. 

The argument above suggests setting up the relation 

2L 

(14) 

so that the matrix Od may be computed explicitly as 

Odr = ^dlrldl). (15) 

Therefore, the duality map Ua is Gaussian if, and only 
if, the matrix Od is invertible, in which case it is also 
orthogonal. In the absence of a more educated choice, 
one may always set 7 ^ = 7 s in Eq. (14). 


The association Ud Od shows that Gaussian duali¬ 
ties, though many-body in nature, also induce a posteri¬ 
ori a duality of the single-particle Hamiltonian h by the 
relation 

2L 

h?s = Y. OdWs,. (16) 

r' ,s' — l 

How much of the locality of the original system one pre¬ 
serves in the dual model will depend on the range of the 
matrix Od- 

Finally, let us contrast Gaussian dualities to other 
type of dualities that do not preserve the quadratic 
fermionic nature of the original theory. Consider the non¬ 
interacting Hamiltonian 

L-l 

H =\'Y% I2jl2j + 1 - ty 72 ^- 172 ^- 1 - 2 ], (17) 

which can be also described as a spin-1/2 Hamiltonian, 

L-l 

77 = - ^ [G CTj crj+i + ty cr^jCr^+i], (18) 

i=i 

after the Jordan-Wigner map of Majorana operators 

721-1=n n ( 1 ®) 

1=1 1=1 

in terms of Pauli matrices ctJ, v = x,y, z. A simple local 
rotation around the spin y axis, 

(j = l,...,L), (20) 

induces a non-trivial change in the dual fermionic Hamil¬ 
tonian. Although local, it is an interacting Hamiltonian 

L-l 

= y][G72j-l72i72i-|-l72j-|-2 - 721 - 1721 - 1 - 2 ] • (21) 

i=i 

Therefore, the matrix Od of Eq. (14) should fail to be 
invertible, and one can check that this is indeed the case. 

The dual Hamiltonian of Eq. (21) has an interesting 
physical interpretation. It describes the competition be¬ 
tween a p-wave superconducting chain in its topological 
phase and a density-density interaction. In the limit in 
which ty vanishes, its ground state is number conserv¬ 
ing, a Mott insulating state, otherwise its ground state is 
superconducting. For G = 0, has two exact zero- 
energy modes, 72 and 72 L-i- The evolution of these 
modes with G may be computed exactly by exploiting 
the duality transformation connecting to the free- 
fermion Hamiltonian H. 


B. Decoupling transformations 

The results of the previous section, especially the ex¬ 
ample at the end of the section, show that generic du¬ 
alities do not preserve the non-interacting character of a 
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theory. Hence, to systematically establish equivalences of 
topological superconductors and insulators it is necessary 
to determine all possible Gaussian dualities. Recall that 
the key difficulty in searching for dualities is to identify 
unitary transformations that respect the locality struc¬ 
ture of the Hamiltonian, meaning that Eq. (14) is not 
the full answer to our problem. We still need to address 
the issue of locality for the specific purpose of relating 
topological insulators to topological superconductors. 

Let us denote by r the sites of a lattice A, i.e., r e A , 
defined in arbitrary space dimensions. A generic (second- 
quantized) electron system where the number of electrons 
N is conserved is described by the Hamiltonian 


H = - 


E 


4 crG 


H” 


( 22 ) 


where c\ ^ represents a canonical fermion creation op¬ 
erator at site r and spin a =t, The hopping am¬ 
plitude tr,r' is related to tr',r by complex conjugation, 
t , = t* , 

^r'.r ^r.r'* 

Generically, the Hamiltonian above displays a broken 
time reversal symmetry unless the hopping amplitudes 
are purely real or purely imaginary. For purely imag¬ 
inary amplitudes, an internal decoupling occurs in the 
system that splits the particle conserving Hamiltonian 
H into four identical, independent and decoupled super¬ 
conductors. The proof of this assertion relies on rewriting 
the particle conserving Hamiltonian H above in terms of 
Majorana fermions ^ and b^. such that 

^r.a ~ G,(T T “ ^r,<T 

with the result 



Hence, if t^y is purely imaginary, H = 
where the particle non-conserving Hamiltonians 


H2,a 


1 

2 

1 

2 


^ ^ tr,r' ^r,(7^r' 
r,r' 

r,r' 


H2,a) 


(25) 

(26) 


of spinless Majorana fermions jr that includes only one 
fourth of the original number of fermionic degrees of free¬ 
dom, and represents any of the four Hamiltonians ob¬ 
tained by decoupling H. 

The equivalence of Eq. (27) states that the spectrum 
of H can be reconstructed from that of i7reduced- In par¬ 
ticular, the zero-energy modes of the particle conserving 
Hamiltonian H are explained by the zero-energy modes of 
the spinless superconductor iJreduced- If I?reduced has pre¬ 
cisely s zero modes, then H has precisely zero modes. 

It is also possible to add spin interaction terms and 
retain some level of decoupling. For example, spin-orbit 
terms of the Rashba, Dimmock, or Dresselhaus type are 
linear in momentum and hence purely imaginary. Thus 
they couple to (and H2,t I® ^ 2 , 1 ), but they 

do not couple a to 5 Majoranas. In the presence of 
these types of spin terms, the decoupling transforma¬ 
tion decomposes the particle conserving system into two 
rotationally-invariant superconductors with non-trivial 
spin dynamics. 

Many of the Majorana lattice models investigated in 
the literature can be interpreted as the ^reduced associ¬ 
ated to a particle-conserving Hamiltonian. The Majo¬ 
rana chain of Kitaev is obtained in one dimension from 
a chain with only nearest-neighbor hoppings.^® In two 
dimensions it is possible to obtain Kitaev’s Honeycomb 
model^^ (see Ref. 28 for the immediate connection) as the 
reduced Hamiltonian associated to the simplest^® model 
of graphene (see below), and the triangular Majorana 
lattice®® from a corresponding particle-conserving model 
on the triangular lattice. Decoupling on the square lat¬ 
tice obtains variations of Majorana arrays investigati- 
gated in Ref. 18. 

An apparently less general but more often useful ver¬ 
sion of the decoupling transformation exits for purely real 
hopping amplitudes on a bipartite lattices A = AU B, 
with generic lattice sites r € A. Let us denote by x G A 
and y G B the sites of each sublattice. The generic bi¬ 
partite Hamiltonian 

H = — ^ ' [lx,y Cx T ^y,x Cy^CT'^x.o-] ) ( 28 ) 

x,y,CT 

allows only for hopping from sublattice A to i? or vicev- 
ersa. If tx,y is purely real-valued, that is, tx,y = iy,x, 
then 


commute, [LIi.o-, = 0. 

The four decoupled superconductors 
Hi^^, H 2 ^ 1 , Hi,^, H 2^1 are isospectral: Any one of them 
can be mapped into any other one by a local unitary 
transformation. On one hand, {H 2 ,-\) is mapped to 
(1^2,i) by a rotation in spin space. On the other 


l + Or abr,, 


hand, the unitary transformation T/o- = Hr ( 

maps Hi a to 772,0-• Hence it is possible to unequivocally 
associate H with a new Hamiltonian iJreduced 


H I y 77reduced — ^ t 


2 tr,r'7r7r', 

r,r' 


(27) 


H = Y, [Hl.a - H2,a] , 

(T 

(29) 

where the superconducting Hamiltonians 


~ ^ ^ ^x,y ^x,cr^y,(Tj 

x,y 

(30) 

^2,(7 2 ^ ^ ^x,y ^XjCr^y,(7 j 

(31) 

commute, [Hi £,, i72,o'] = 0. The spectral equivalence of 
Hi a- and 772 ,0 is established by the unitary transforma- 
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tion 

7 y _ I I f ^ ^X,(T^X,(T \ f ^ \ 

""-iil V 2 ) I V2 ) 


(32) 


that maps i/i^cr ^ -ff 2 ,CT- An interesting corollary to Eqs. 
(29) and (32) is that the unitary transformation C = 
hl-\hli anticommutes with the Hamiltonian, CH = —HC, 
and so it defines a chiral symmetry in the sense that the 
spectrum of H is particle-hole symmetric, i.e., for each 
positive eigenvalue Ea there exists a negative —Ea- 

Just as before it is possible to associate a reduced 
Hamiltonian to El. 

Consider, as an example, the chain of spinless fermions 


2M 


H = 


■ ^ tj - 1/2) . (33) 


i=i 


with quenched disorder in the real-valued hopping am¬ 
plitudes tj and on-site atomic energy tj. The chain has 
L = 2M lattice sites and periodic boundary conditions 
are assumed = c|). Since the lattice is bipartite, 

H is the difference of two identical, independent super¬ 
conductors, coupled by the on-site atomic energies. Let 
us associate a pair of Majorana fermions aj,bj to each 
site j, as in Eq. (23), and rewrite the Hamiltonian above 
in terms of Majorana degrees of freedom 

. 2M 

H = --J2[t3 (aj+ibj - bj+iaj) + e, a^b^] 

i=i 

= ill - + H,, (34) 

with commuting Hamiltonians 
i ^ 

Hi = “2 X! a2jb2j-i - t2j b2j+ia2j] , 

i=i 

. M 

iJ2 = “2 X! b2ja2j-i - t2j a2j+ib2j ], (35) 

i=i 

and 

M 

He = -(i/2) ^[e 2 j-ia 2 y-i 62 j-i -f e2ja2jb2j]- (36) 
i=i 


C. A class of Gaussian duality transformations 

We are now ready to introduce a large class of Gaussian 
duality transformations. Consider the cases for which the 
sublattices A and B of previous section are equivalent, 
meaning that there is a shortest, typically non-unique 
translation (5i (^ 2 ) mapping sublattice A (B) to sublat¬ 
tice B (A). In set notation, 

A + Si = B, B + d 2 =A. (37) 

For the purpose of the duality transformation that we are 
about to introduce it is often convenient to choose ^ 1 ,62 


to be as parallel and short as possible. This condition 
guarantees that the range of the hoppings in the dual 
Hamiltonian deviates as little as possible from that of the 
original Hamiltonian. A hypercubic lattice is simplest in 
that one may choose = ^ 2 - 
The mapping 


^y,(7 ^ ^X,(7 ^ ^X.CTl 

^y,(7 ^ by fjj b^ fj y (38) 

induces a unitary transformation that leaves Hi ,j un¬ 
changed and transforms 772,0- as 

772,0 ^ -7^2,0 — 2 ^ ^ ^y—5i,x—(52 7x,oay,o. (39) 

x,y 


In rearranging the sum over sites, we have assumed pe¬ 
riodic boundary conditions or that the system is infinite. 
The dual superconducting Hamiltonian, 

77"" = E [77i,o - 77/^.] = (40) 

<7 

X ^ ^ [^3C,y ax,o7y,o 7y_5y^X—62 ^X,Oay,o] ( 
x,y,o 


can be rewritten in terms of creation and annihilation 
operators. 




- E [^^^(ci.oCy.o + 4.-Cx,o) + 

x,y,^T 

Ax,y(c)y^^c)(.^^ -I- Cx^erCy ,y)j , 


(41) 


where 


j.av _ 

^x,y — 


7 x,y 4 “ 7 y_^^ X —62 


A _ ^x-y ^y- 

^X,y - 2 


5i,x—(S 2 


Even though 77 J- 77^ is roughly as general as possible 
for a band electronic system of independent fermions, 
77^ + remains a superconductor at vanishing chem¬ 
ical potential fx. It is possible to include spin terms in 77 
and still obtain a dual superconductor featuring only lo¬ 
cal interactions. Just as the duality breaks particle con¬ 
servation in general, we expect it to modify rotational 
properties since it has a highly non-trivial action on the 
operators of total spin. Since, and when, spin does not 
play any decisive role in the studied physical phenomenon 
we will drop it from the discussion in order to avoid con¬ 
fusing notation and obscure explanations. 

It is now straightforward to apply the general dual¬ 
ity transformation of bipartite models to the disordered 
chain of the previous section. For this one-dimensional 
system, the mapping defined in Eqs. (38) reduces to 


o-2j-i —> b2j, a2j —>■ a2j 

b2j-i —> b2j-i, b2j —>■ —a2j+i, (42) 

always identifying the index L + 1 with I, and 0 with 
L — I. Thus, while Hi = TJf* remains invariant, H 2 
transforms as 

. M 

H 2 = “2 E [^2i-2 b2ja2j-i - t2j-i a2j+ib2j] (43) 
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(the case M = 1 is special in that the duality map keeps 
H 2 also invariant). The on-site atomic energy term 
transforms like 

. M 

He = “2 ^ [^ 2 j-lh 2 jh 2 j-l - £ 2 ^ 02 ^ 02 ^ + 1 ] (44) 

. 2M 

= ^ H - c]+iS) + - c^-+ic^) . 

i=i 

Combining all of these results, we obtain the dual super¬ 
conductor = Hi- 

1. Symmetry transmutation: particle number and fermionic 

parity 

The duality transformation Eqs. (38), breaks parti¬ 
cle conservation in general because the particle number 
(charge) operator 

H = '^[nr,c, -^/‘A = ( 46 ) 

r.C r,(T 

associated to, and a symmetry of, H is drastically mod- 
ihed by the duality. Since 

Ox,(T^x,cr t Ox,ctOx+5i ,cr 5 

t^y,aby^(7 t ^y,o'^y+52,f7; (4^) 

the duality transformation maps iV to a symmetry 
of that does not have the interpretation of a charge 
operator, 

N = —— ^ ^ ] Qx,a-Qx+(5i,o- + ^ ' ^y,(T^y+^2,o’ j (48) 

^ L X y 

while it is still true that [N^, H^] = 0. 

There is, however, a quantum number important from 
the point of view of superconductivity that is almost pre¬ 
served by duality: fermionic parity. The operator of 
fermionic parity 

(-1)"’ = = n(-iar..6r..) (49) 

r,(T 

measures the parity of the total number of fermions. The 
BCS mean field approximation breaks the symmetry of 
particle conservation down to conservation of fermionic 
parity. The duality transformation maps 

(-1)^ ^ (50) 

J_(io.x,crO-x+6i,(T) J^(i5y,(T5y+52,o-) = (~1) (~1) ) 

x,(T y,a 


where (—l)'^ is the sign accumulated after permutations 
of the Majorana fermions to establish the original order 
(- 1 )^. 

Incidentally, Eq. (47) shows that the duality H —> 
is more general in scope than the decoupling transforma¬ 
tion that motivated it. For example, adding an on-site 
energy term 

He = -J2^rinr,a-l/2) (51) 

r.cr 

(r G A U B) to H, couples the reduced superconductors 
Hi er and 772,o-- It transforms as 

He = 2 ^ ( £xQx,iTQx+3i,cr + “ ^ ' £y5y.cr^y+ga,o-■ (52) 

x,cr y,a 

Hence the effect of the on-site atomic energy term Cr is to 
renormalize (by a purely imaginary amount) the hopping 
and pairing amplitudes of the dual superconductor in the 
directions (5i, and ^ 2 - 

2. Translation Symmetry 

For the Gaussian dualities of Eq. (38), the non¬ 
conservation of particle number for the dual partner 
is precisely related to (partial) breaking of translation 
symmetry for H, since the pairing potential Ax,y van¬ 
ishes if ty-Si,x-S 2 = 7x,y (and = 77 in this case). 
What is less obvious is that the translation symmetry of 
77^ may be be higher than that of 77, in which case we 
are enlarging one group of symmetries (translations) at 
the expense of breaking another symmetry, particle con¬ 
servation. This is explained by the transmutation of the 
translation operation under duality. 

Let us focus for simplicity on a closed chain of spin¬ 
less fermions. The extension to more general settings is 
straightforward but notationally cumbersome. The map 
of Eq. (42) leads to 

C2j-1 c^_i = + \h 2 j-l), 

C 2 j ^ eg = ^{a 2 j - \a 2 j+i), (53) 

or, more explicitly, 

cg-i = ^(C2i-i - ic2i) - \ - i4i): 

eg = ^(C2i - ic2j+i) + ^(cg - icg+i). (54) 

These expressions show already transmutation of the 
translation operation, let us call it T. On one hand, 

fc,ft=c,+i {j=j + L), (55) 

and consequently {T^ =U^TU\), 
f^cg(f^)t = cg+i 


{j = j + L)- 


( 56 ) 



On the other hand, 


3. Time-reversal Symmetry 


(57) 

and so 

f^f^. (58) 

This is the point to notice. If T happens to be a sym¬ 
metry of H, then is necessarily a symmetry of H^. 
However, cannot possibly have the interpretation of a 
translation by one site, since that physical interpretation 
continues to be attached to T ! (The action of on 

the Cj may be computed by inverting Eqs. (54).) Notice, 
however, by the same reasoning that 

^ giad(fr')2^ (59) 

where the (possibly trivial) phase on the right-hand side 
is determined by the actual duality transformation. 

The concrete signihcance of this result will become ap¬ 
parent in the next section when we investigate the equiv¬ 
alence of the dimerized Peierls chain and the Majorana 
chain of Kitaev. What happens in that case actually is 
the following. The dimerized Peirls chain commutes with 
and a very non-evident symmetry u\ T (not to be 
confused with =U^TU\). Asa consequence, its dual 
partner (the Majorana chain) commutes with T. This il¬ 
lustrates how Gaussian dualities may increase translation 
symmetry at the expense of breaking (transmuting) other 
symmetries, e.g., particle conservation. 

It is revealing to to rewrite the Gaussian duality of Eq. 
(42) as a map of fermions in crystal momentum space. 
The even-odd structure of the duality mapping evinced 
by Eqs. (54) for example suggests that we should take a 
unit cell with two sites. To keep the notation simple, we 
will assume that L = 2M, with M odd. Then 

M-1 M-1 

C2j-i= ^ ^ ~ 7 W 

M-1 M-1 V 

^— 2 ^— 2 

where k = 2 ttI/M. Let us emphasize that we are not 
making any assumption about the symmetries of any par¬ 
ticular Hamiltonian. We are just going to recast our dual¬ 
ity transformation in a new light. With these definitions, 
the dual fermions in momentum space are 

Cpfc = ^(ci,fc - ico.fc) - (61) 

= ^(co.fc - \e-'^ci,k) + ^(4.-/0 - 

The key point is that the dual fermions of momentum k 
are combinations of the original fermions of momentum k 
and —k. It follows that the induced single-particle dual¬ 
ity Od is not block diagonal with respect to momentum. 


The standard antinnitary operation T of motion re¬ 
versal may be specihed by its action on the creation and 
annihilation operators, 

'1' S-.T = df’ 

'1' S-4 4; = -4t> (62) 

with the important consequence that 

= (-1)^. (63) 

The duality transformation Eq. (38) is spin diagonal. 
Hence the dual fermions 

C2j-1,cr = 2 ('"21-1.'^ ~ '^ 2 ^, 0 -) — ~ 

^2j,<T — 2^C2j^a — iC2j-|-l,CT) + 2^C2j,<7 ~ + (64) 

are just as before, except for the additional spin label. 
By construction, the dual antiunitary operation = 

UdTU\, with 

(r^)2 = (_i)^", (65) 

acts as standard time-reversal on the dual fermions (for 

a discussion of , see Section HCl). One may check 
that 

r^r^, (66) 

and so there is transmutation of time-reversal symmetry. 
To put this result in perspective, suppose that both H 
and commute with T. As we will see, this is the case 
for example in polyacetylene and its dual superconduct¬ 
ing partner (class Dili). Then, since [i7^,T'°] = 0, we 
have uncovered a (unitary) symmetry T^T of . 

An example of transmutation of time-reversal for 
canonical bosons (phonons) can be found in Ref. 11, 
page 730. 

III. EQUIVALENCES OF TOPOLOGICAL 
INSULATORS AND SUPERCONDUCTORS 

Can Gaussian dualities in general, and in particular, 
the ones of this paper, establish equivalences between 
topological insulators and topological superconductors? 
The discussion of the previous section shows that there 
is no obstruction for this to be the case: many-body du¬ 
alities can jump across entries in single-particle classifi¬ 
cation schemes simply by transmuting key symmetries. 
But symmetry transmutation is a necessary, not a suf¬ 
ficient condition. For example, the duality of Eq. (42) 
maps the clearly trivial insulator, 

2M 

H = -e^(n,-l/2) 


(67) 
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to the equally trivial superconductor, 

. 2M 

= f H [(S^S+1 - c]+iS) (68) 

i=i 

+(-l)^(c]c]+i-Cj+iCj) , 

in spite of the symmetry rearrangements it causes. 

For the Gaussian dualities of the previous section in 
particular, it is not hard to convince oneself that it must 
be the case that they map topologically (non)trivial sys¬ 
tems to equally non(trivial) dual partners. In this section 
we will study some paradigmatic examples of non-trivial 
partners. In one dimension, we find that the dimerized 
Peirls chain and the Majorana chain of Kitaev are dual 
partners, and we also investigate the mapping of topolog¬ 
ical defects under duality. The more general m—merized 
Peierls chain is investigated in the Appendix A. In two 
dimensions, we study a topological insulator based on a 
Kekule-like pattern of hopping matrix elements that in¬ 
cludes graphene as a special case. It dual partner is a 
p-wave superconductor. In the limit where the insulator 
becomes graphene (a semi-metal), the dual superconduc¬ 
tor reduces to a stack of Kitaev chains interconnected by 
pure kinetic hopping in the direction perpendicular to the 
chains. This superconducting realization of Dirac cones 
seems to be new in the literature. Appendix B shows an 
equivalence of a trivial BCS superconductor to a trivial 
insulator in any number of space dimensions. 


A. The Peierls chain is dual to the Kitaev chain 

In one dimension, the dimerized Peierls chain at half¬ 
filling, proposed by Su, Schrieffer and Heeger (SSH) to 
model polyacetylene^^, is the prototype of a topologi¬ 
cally non-trivial band insulator, while the Kitaev chain 
is the prototype of a topologically non-trivial supercon¬ 
ductor. In spite of their physical differences, the Kitaev 
and Peierls chains are isospectral in second quantization, 
that is, as many-fermion systems. The reason is that 
there exists a Gaussian duality connecting both models. 
As we will see below, the mapping is different in detail for 
periodic and open boundary conditions. The Gaussian 
duality for open boundary conditions is crucial to under¬ 
stand the way boundary excitations are related, i.e., how 
the Majorana charge-neutral zero-energy edge modes of 
the Kitaev chain map into charged (canonical fermion) 
zero-modes of the Peierls chain. 

For simplicity, and pedagogical reasons, in the follow¬ 
ing we will consider the spinless case. (The original model 
for polyacetylene^® involves spin-1/2 electrons and this 
fact is relevant for topological classification purposes but 
it is not from the standpoint of Gaussian dualities that 
preserve spin, such as the ones defined in this paper.) 
Since the unit cell of the dimerized Peierls chain consists 
of two sites, the length (number of sites) L of the closed 
chain must be even, L = 2M. There are only two in¬ 
dependent, periodically repeated hopping terms ti and 


j = I j = 1 

L = 2M t2 h 2 L = 2M ti2 




FIG. 1. Peierls chain with periodic boundary conditions on 
the left and its dual Kitaev chain superconductor on the right. 
Here, L = 2M. 


^ 2 , see Fig. 1. It follows that the Peierls chain is a spe¬ 
cial case of the generic one-dimensional Hamiltonian of 
Eq. (33), with 


t2 j-l=ti, t2j=t2, €j = 0, j = 1,-■ ■ ,L. (69) 
The dual of the Peirls chain, 

L 




o.j+ibj — t2 bj+ic 


(70) 


i=i 


L, 

■n 

i=i 


ti2 c]c^+i + Ai 2 , -I- H.c. 


DD+1 


is obtained by the same specialization of Eq. (45), with 

f . _ti-t2 , . 

tl2 — -^-, Ai2 - - - -. (7ij 

The Hamiltonian is precisely the Majorana chain of 
Kitaev, at vanishing chemical potential fa. The interplay 
of symmetries is noteworthy. While the dimerized Peierls 
chain shows reduced translation symmetry and conserva¬ 
tion particle number, its dual the Kitaev chain has full 
translation symmetry and particle-number conservation 
is broken. 

Next, we would like to comment on polyacetylene 
(class AH), that is the case with real electrons (spinfull 
fermions). In this case, we obtain two copies of a Peierls 
chain, one for each spin component. Time-reversal sym¬ 
metry is preserved and our Gaussian duality maps the 
SSH model to two copies of Kitaev’s chain (class Dili), 
one for each spin component. 

In Appendix A we analyze the superconducting equiv¬ 
alent of an TO-merized Peierls chain with m > 3. There 
is a clear distinction between distortions with m even 
and odd. In the latter case, the periodicity of the dual 
superconductor is doubled, and there are no zero-energy 
modes. For m > 4 even, the periodicities of the insulator 
and its dual superconductor are the same, in contrast to 
the dimerized case discussed above. It is well-known that 
for m even there are zero-energy modes. 
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1. Mapping of topological invariants 

A main goal of any topological classification of matter 
is to divide all possible quantum states of matter into 
equivalence classes. There is some degree of arbitrari¬ 
ness in the criteria used to define those classes. Once the 
criteria is established, e.g., by symmetry/dimension and 
reduced K-homology of bundles,^ two states in the same 
class are connected through a continuous map whose in¬ 
verse is also continuous, i.e., a homeomorphism. To es¬ 
tablish a characterization of the class one uses topological 
invariants, i.e., quantities that are preserved under the 
homeomorphism. That a particle-conserving system such 
as a Peierls insulator may be dual to a superconductor 
raises some conceptual issues for the topological classifi¬ 
cation of systems of free fermions. What is the relation 
between the topological invariants characterizing these 
different but dual states of matter? The fermionic parity 
of the ground state of a superconductor such as Kitaev’s, 
with L = 2M even, constitutes a good quantum num¬ 
ber whose value depends on the boundary conditions. 
For periodic boundary conditions, the topologically non¬ 
trivial ground state of the Kitaev chain is non-degenerate 
and fermionic parity is odd, while it is even in the triv¬ 
ial phase.These facts hold independently of M = T/2 
(the relevance of this remark will become clear below). 

Fermionic parity may be computed in several equiv¬ 
alent ways. The topological character of this quantum 
number is revealed by its connection to the Majorana 
number,^® 

= sgn(^-k2ti2) sgn(/r- 2 ^ 12 ), (72) 

defined as the product of signs of Pfaffians of an anti¬ 
symmetric matrix at momenta 0 and -tt. This quantity 
identifies the topologically non-trivial phase as the one 
with = — 1. On the other hand, we have seen that 
the dimerized Peierls chain maps into a Kitaev’s chain at 
p. = 0. The Peierls chain is always in a topologically non¬ 
trivial insulating phase, as long as ti 7 ^ O- However, the 
fermion parity of the Peierls’ insulating non-degenerate 
many-body ground state is given by (—1)^, and therefore 
it is defined by the parity of M, i.e., it can be odd or even 
depending on M. It is instructive to express the Peierls’ 
chain in terms of Majorana fermions and compute the 
Majorana number, now with a doubled unit cell, to real¬ 
ize that indeed A4 = (—1)'^. The point is that fermion 
parity is not the good topological quantum number to 
characterize the Peierls’ insulating phase despite the fact 
that it is exactly dual to Kitaev’s chain model. 

From the point of view of the many-body duality trans¬ 
formation, the mismatch is explained by the (very mild) 
transmutation of fermionic parity, Eq. (50). 

2. Mapping of topological defects and boundary modes 

Consider for simplicity a dimerized, spinless, chain 
with periodic boundary conditions, L = 2M with M G 


odd, and two defects symmetrically located at positions 
j = 1 and j = M. This corresponds to the Hamiltonian 

M-l 

H = — ^ —l(mod 2 ) (CjCj_|_]^ (73) 

L-1 

— ^ ij-(M-|-l)(mod2)(c]Cj + i + c]_,_;^Cj-) 

j = M+2 

— t2{c\c2 + + C^M+1^M+2 + H.C.), 

where the first two terms represent two identical dimer¬ 
ized Peierls chains each of length M — 1, and the last 
term represents the pair of defects, see Fig. 2. 




FIG. 2. Peierls chain (periodic boundary conditions) with a 
couple of defects on the left and its dual Josephson junctions 
of Kitaev chain superconductors on the right. Here, L = 2M. 

Applied to the Hamiltonian (73), the duality trans¬ 
formation of the periodic chain, Eq. (42), produces two 
dual superconductors coupled by two particle-conserving 
segments, 

M-l 

= -Yl + H-C-) (74) 

L-1 

- (^I2cjc^ + 1 - Ai2c];c] + i -h H.c.) 

j^M+2 

-t2{c\c2+C^M+l^M+2+^-C-)- (75) 

Notice the change in sign of the superconducting order 
parameter across the links. Because of this phase dif¬ 
ference, a Majorana zero mode is trapped at each weak 
link. 

For an infinite chain with a single defect located at 
the origin, one may use the same duality map, Eq. (42) 
in order to obtain a dual Majorana zero-mode localized 
at the origin. The defect that famously traps fractional 
charge ±e/2, per spin direction, in the Peierls chain is 
dual to a defect that traps a Majorana zero-mode! 

For open boundary conditions, it is necessary to take 
the number of sites L = 2M -|- 1 to be odd, see Fig. 3. 
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The duality transformation 

a2j-i —>■ ai_2(j_i) , a2j —>■ a2j 

b2j-i —>■ b2j-i , b2j —>■ (76) 

leaves Hi = invariant, and transforms H 2 as follows 

i ^ 

H 2 = “2 aL-2jbL+l-2j - t2j-l &L+l-2jaL-20-l)), 

i 

= ~2 ®2j-1^2j ~ ^L-2j ^2ia2i + l)- (77) 

^ 1=1 


Let rj = t 2 /ti, and assume without loss of generality |? 7 | < 
1. From these commutators it is possible to show that 
the combination 

1 

^ “2j+l (81) 

1=0 


of Gj fermions at odd sites j commutes with the Hamilto¬ 
nian . It constitutes an exact symmetry for any finite 
M. The normalization factor is 


Af{r]) 


1 _ ^2(M+l) 


1 _ 


(82) 


tl t2 t\ ti t2 

•=• - • - • ■ ■ • - • - •=• - • 

j = 1 2 • • • L = 2M + l 

tl2 ti2 ti2 ti2 ti2 

• - • - • - • ■ ■ • - • - • - • - • 

Ai2 Ai2 Ai2 Ai2 Ai2 

FIG. 3. Peierls chain with open boundary conditions on the 
top and its dual Kitaev chain superconductor on the bottom. 
Here, L = 2M + 1. 

In the dimerized case there are only two different al¬ 
ternating hopping terms which satisfy 

t2j-l = tL-2j = tl , t2j = tL+l_2j = t2, (78) 


A similar calculation establishes the right symmetry 

1 ^ 

7r^ht = JfYj - (83) 

1=0 

These exact Majorana zero-energy modes are exponen¬ 
tially localized. If \ri\ > 1, the corresponding localized 
symmetries are obtained by rescaling 7 ^ —1 
and changing the normalization factor to JV{l/ri). At 
\rj\ = 1 the mass gap vanishes in the thermodynamic 
limit L —1 00 . 

Notice the fundamental difference between chains of 
even or odd lengths L. While it is possible to determine 
exact zero modes when L is odd, this is not the case for 
L even where the Majorana character of the edge modes 
is only asymptotically exact in the thermodynamic limit. 
The reason is simple. For any finite L, commutes 
with the global symmetries 


with the end result that the dual total Hamiltonian rep¬ 
resents a spinless superconductor 

. L-l 

= -^ ^(H aj+ibj - t2 bj+iOj) (79) 

1=1 
L-l 

= - ^ (ti 2 + Ai 2 + H.c.) , 

1=1 

that is again the Kitaev chain Hamiltonian at vanishing, 
/i = 0, chemical potential. 

The Kitaev chain of this section (open boundary condi¬ 
tions, odd length, and vanishing chemical potential) has 
two exact zero-energy modes, one per boundary point. 
The chain is reflection symmetric with respect to the cen¬ 
tral site j = M + 1, and the many-body ground state is 
two-fold degenerate. The zero-energy mode associated 
to the left boundary can be computed from the set of 
commutators 

[-\H^,ai] = t 2 b 2 , 

[- 177 ^, 03 ] = ti&2 + ^2^4, 

[—177^, aL_2] = tihr-d, + 72&l-i, 

[-177^, ol] = tibr-i- 


L L 

U, = 1[{1- 2n,) , [/, = n Kb] + &,), (84) 

1=1 1=1 

where b] = c] (1 ~ ^nj) is a hard-core boson. How¬ 
ever, it is only when L is odd that {Uz,Ux} = 0. In 
turn, this implies that the whole many-body spectrum of 
77^ is, at least, exactly two-fold degenerate. This sym¬ 
metry analysis applies just as well to the more general 
Hamiltonian of Eq. (33) with open boundary conditions, 
provided the on-site potential Cj vanishes. 

The duality transformation of Eq. (76) maps the 
boundary Majorana zero-modes of the Kitaev chain into 
corresponding boundary symmetries of the Peierls chain, 

1 “ 

7right,l = 'y ( (~ll) 0'L-2j, 

1=0 

1 “ 

7right,2 = JfYl ^L- 2 i- (85) 

1=0 

Unlike for the Kitaev chain, these two zero modes reside 
on one and the same edge. Hence, it is natural to re¬ 
combine them into one exponentially localized fermionic 
mode, 

t _ 7right,l — i7right,2 
right 


(80) 


2 


( 86 ) 













12 


It is not surprising that both boundary symmetries of the 
Peierls chain appear on one boundary point, and not both 
as in the Kitaev case. The reason is the lack of reflection 
symmetry about j = M + 1 in the mapped Peierls chain. 
Nonetheless, its many-body ground state is two-fold de¬ 
generate, just as its dual Kitaev superconductor, indicat¬ 
ing that the Z 2 symmetry of fermionic parity is odd in 
the thermodynamic limit also for the Peierls chain. This 
statement is of course confirmed by the exact solution of 
the Peierls Hamiltonian. 


3. Density-density interactions 


As mentioned above Gaussian dualities can be used to 
establish equivalences in interacting many-body systems. 
Here, we illustrate this fact in another paradigmatic ex¬ 
ample. Consider the case of a dimerized Peierls chain (of 
length L = 2M) at half-filling where electrons interact 
through a density-density Coulomb repulsion V 


2M 


H = 




+ 4 + 1 C) 




- C- [CjC^ + l -r 


-|- 4 (1 2nj){l 2nj+i) 


(87) 


The duality transformation of Eq. (42) maps the density 
operators as follows 


n2j-i 2 ^^ ~ 'hj-ihj), 

n 2 j -t 2(1 “ 102^02^+1), ( 88 ) 

and the resulting dual superconducting equivalent is 
given by 

L 

77^ = - ^ [<12 c^Cj+i + Ai 2 -t H.c. 

V ^ 

-t— V 
4 

J=i 

which clearly shows the competition and interplay be¬ 
tween the band and Mott gaps. 

A natural question that emerges is how robust is the 
Peierls phase to the presence of Coulomb interactions? 
The second, related, question is can interactions alone 
generate a topological Mott phase in the case where the 
non-interacting phase is metallic, i.e., A 12 = 0? Since 
in the latter case the model is exactly (Bethe ansatz) 
solvable, we know that for sufficiently large repulsion V, 
there exists a Mott phase. 

In Appendix C we analyze the phenomenon of symme¬ 
try transmutation in interacting boson systems. 


„t„ I Jrj t„t 

SS-^2 + e C^C^+2 


(89) 



ex 





t'l 


t '2 




h 

t2 



X 


y 








Dual Superconductor 





r' 




j.av 

tj, j./ 

zXj. j./ 




j.av 

r 


— 77 

r 

+77 

zXj, j,/ 

r' 


+77 

— 77 





FIG. 4. (Top) Topological insulator characterized by a par¬ 
ticular Kekule-type pattern of hopping matrix elements as 
depicted in the figure. Double bonds along the horizontal 
(vertical) direction represent the hopping matrix element ti 
(t'l), while single bonds along the horizontal (vertical) direc¬ 
tion represent t 2 (t^). The ex{y) are on-site atomic energies, 
constant on each sublattice. (Bottom) Dual chiral topolog¬ 
ical superconductor with tt fluxes per square plaquette dis¬ 
tributed antiferromagnetically, and a -|- ipy superconduct¬ 
ing order parameter. The dual relation between parameters 
is explained in the main text. 


B. Graphene is dual to a “weak” Topological 
Superconductor 


In this section we establish a non-trivial equivalences in 
two space dimensions. For conciseness, we investigate a 
topological insulator on a square lattice, characterized by 
a Kekule-like pattern of hopping parameters 
and an on-site potential ex,ey that is constant on each 
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sublattice, see Figure 4. In this way we manage to ad¬ 
dress several interesting models in a unified fashion. The 
semimetal graphene for example is realized on the line 
ti = =< 2)^2 = 0 = Cx- (A few other lines obtain 

graphene as well. Notice that the honeycomb lattice is 
appears represented as a brick wall lattice in Figure 4). 
As a consequence,^^ our model realizes a condensed mat¬ 
ter analog of the (2-|-l)-dimensional parity anomaly^^, 
and our duality transformation provides a superconduct¬ 
ing dual representation of this phenomenon. As was the 
case in one dimension, the Gaussian duality mapping the 
insulating model to a topological p-wave superconductor 
has the key property of preserving the locality of the edge 
mode excitations. 

The two-dimensional underlying lattice A considered in 
the following is bipartite with lattice points x = (xi,X2) 
and y = {yi,y 2 ), such that Xi + X 2 € even and yi + 
y 2 € odd integers. The total number of lattice points 
along the horizontal direction is La,, and Ly along the 
vertical direction, such that Lx x Ly defines the size of 
the lattice and where, for simplicity. Lx and Ly represent 
even integers. Figure 4 (Top) is an example of a lattice A. 
Consider in particular the lattice shown in Fig. 4 (Top) 
where, for a given point x = (xi,X2), the corresponding 
hopping amplitudes of Hamiltonian (28) are given by: 

{ h , for y = (xi + 1 ,X 2 ) 
t 2 , for y = {xi - 1 ,X 2 ) 
t'l , fory=(xi,X 2 -l) 
t '2 , for y = {xi,X 2 -f 1) 

and with on-site energies Cx = —Cy (see Eq. (51)). 

This model, endowed with periodic (toroidal) bound¬ 
ary conditions, has a single-particle energy spectrum 
(bulk bands) given by 



/^x + + + Bl _ , 

1 -^4,k(T — -^l,k(T5 

-^2,k(T — 

y^x + Ak _ + : 

1 As^ko- = —E2,kcr,(9I) 


where the wavevectors k = {kx,ky) are defined in the 
Brillouin zone {kx = ^rix, ky = ^riy) with Ux = 

0, 1 , • • • , ^ — 1 and riy = 0, 1 , • • • , ^ — 1, and 

Ak,± = (H + t 2 ) cos ± (t'l -b t' 2 ) cos ^y^ 

^k,± = (H - t 2 ) sin ^y ^ ± (t) - t' 2 ) sin .(92) 

There is a chiral symmetry at work, since the energy 
levels are symmetrically distributed around zero energy 
and time reversal is not broken. 

The Gaussian duality of Eq. (38) with = <^2 = (0,1) 
maps our Kekule-type insulator into a p-wave supercon¬ 
ductor with an antiferromagnetic distribution of 7r-fluxes 
per square plaquette (or more precisely, two spin copies of 
this system). The dual (chiral) superconductor is shown 
in Fig. 4 (Bottom) and corresponds to (see Eqs. (41) 


and (52)) 


H^ = - C'4.^G'.^ + (C')*4'.^G.^ + 

(r,r'),o- 

c]., + (Aj-y) Cr^o-Crgo-)] ) (93) 


where (r, r') represents nearest-neighbor links of a rect¬ 
angular lattice with lattice points r = (ri,r 2 ) and 

for r' = (ri -|- l,r 2 ) 
for r' = (ri,r 2 -b 1), 


j.av _ 

6r r' — 


ti -\-tL 


tl-\-t2 

2 

i(_l)'-l+^2 


A 


r.r' 


ti—t2 



for r' = (n -b l,r 2 ) 
for r' = (ri,r 2 -b 1) 


.(94) 


1. Mapping of topological boundary modes 

Our Kekule-type insulator may displays zero energy 
modes if the on-site potential vanishes. It is instruc¬ 
tive to consider explicitly the case of graphene to high¬ 
light the differences between zig-zag and armchair edge 
terminations^^ from the point of view of the dual super¬ 
conductor. So let us take open boundary conditions along 
the ri-direction and periodic along the r 2 -direction, i.e., 
and open cylinder, and the parameter set G = t) = t' 2 , 
^2 = 0 = Cx = Cy. Then, the Gaussian duality map 
used above for the toroidal boundary conditions (bulk) 
also works for the cylinder, since = ^2 = (0,1) de¬ 
scribes a translation along the periodic direction. One 
can see from Fig. 4 that this situation corresponds to 
a zig-zag edge, while the parameter set G = <2 = t'l, 
t '2 = 0 = Cx = £y would correspond to an armchair ter¬ 
mination. The corresponding dual superconductors rep¬ 
resent (two copies of) a stack of horizontal or vertical Ki- 
taev chains respectively, in the topologically non-trivial 
regime. When the chains are horizontal, they obtain 
the topological superconducting Majorana edge modes 
dual to zig-zag terminated graphene. When the chains 
are vertical, the superconductor does not display edge 
modes, and its dual corresponds to the armchair termi¬ 
nated graphene. 

The general duality transformation of Section IIC 
breaks down for open boundary conditions in both ri and 
r 2 directions. Nonetheless, all of our conclusions hold just 
as well in this case were zig-zag and armchair termina¬ 
tions coexist. In order to illustrate this point explicitly it 
becomes necessary to introduce a different Gaussian du¬ 
ality, showcasing once more the fact that exact dualities 
are very sensitive to boundary conditions. 

Let us focus for simplicity on a particular, spinless, case 
of our Kekule-type insulator with open boundary condi¬ 
tions in both directions. The Hamiltonian of interest is 
given by 


H = H 1 + H 2 + HX, 


(95) 
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where 


Hi 


My 


t 

2 ri —l,2r2 —1 


r2 — l ri=l 


C2ri,2r2-1 


+^2 C 


t 

2ri,2r2 —1 


‘'2ri + l,2r2-l 


H.c. 


(96) 


H2 


My 


Y. Y [^2(4n_i. 


-I,2r2 ^2ri,2r2 


r 2 — l ri—1 


+^1 C2ri,2r2''2ri + l,2r2 + 


(97) 


[i7i,i72] = 0, and 


Ly 1 Lx 

H.=-tlYYblur.Cy^ ,r 2 + l + H-C.j , 

r 2 — l ri—1 


(98) 


with Lx = 2 Mx + 1 and Ly = 2 My, has an appealing 
interpretation as a stack of dimerized Peierls chains. Be¬ 
cause the Peirls chain are alternating, the bulk transla¬ 
tion symmetry in either direction is generated by trans¬ 
lations by two sites. As a function of t 2 , the stack in¬ 
terpolates between a trivial metal, t 2 = ti, and spinless 
graphene, t 2 = 0, with zig-zag vertical and armchair hor¬ 
izontal boundaries. 

The description of as a stack of Peierls chains sug¬ 
gests a natural way to map the system to a superconduc¬ 
tor. Let us apply, to each horizontal chain, the duality of 
Section III A 2, Eq. (76). Unfortunately, this simplest al¬ 
ternative would obtain a non-local dual representation of 
Hy. There is, however, a way to fix this problem. Recall 
that the Gaussian dualities of this paper are motivated 
by the observation that some systems may be split into 
independent subsystems, and then it is possible to rear¬ 
range one subsystem relative to the other. So, given the 
splitting of the open Peierls chain into two subsystems, 
we could rearrange one subsystem for, say, the chains 
at odd height 2 r 2 — 1, and the other subsystem for the 
chains at even height 2r2. This idea is implemented by 
the Gaussian duality 


a2ri-1.2r2-l _2(ri-1),2r2-1) 

^2ri,2r2-l ^Lj,-|-l-2ri,2r2-1) 

^2ri,2r2 ^ —2ri,2r2 5 

^2ri-l,2r2 6L,„-2(ri-l),2r2 • (99) 

The Majorana operators that are not explicitly listed re¬ 
main unchanged. 

Now, the alternating structure of the duality mapping 
in the vertical direction obtains the trivial transformation 

H^=Hy. (100) 

The effect of the first half of the transformation on Hi 
follows immediately from the work in Section III A 2. In 
H 2 , the hoppings ti and ^2 are exchanged relative to Hi. 
However, the second half of the transformation is also 


modified relative to the first half in such a way as to 
precisely compensate for this exchange. Explicitly, 

H^ = H^ + H^ + Hy, (101) 


with 


iff + H^ = 


'Y^ 'Y^ |^G2c4,r2'^'’l+l.'r2 

i’2 = l ri = l 


+ ^124i.r2 4i-|-l,r2 + H-C. .(102) 


As before, the dual system may be described as a stack 
of Kitaev wires, but now with open boundary condi¬ 
tions in both directions. It realizes on its vertical bound¬ 
aries the Kitaev edge,^® a one-dimensional p-wave super¬ 
conductor robust against statistical translation invariant 
disorder and/or interactions.^® We see that the zig-zag 
boundary of spinless graphene is dual to the Kitaev edge. 


IV. SUMMARY AND OUTLOOK 

In this paper we have developed the general theory of 
Gaussian duality transformations for fermions, defined 
as maps that preserve a quadratic Hermitian form. The 
dual partner of a system of free fermions is also a system 
of free fermions if the duality transformation is Gaussian, 
and both systems are equally local in space. The the¬ 
ory and practice of fermionic Gaussian dualities benefits 
from “Majorana fermions,” the complex Clifford algebra 
canonically represented in Eock space. As a consequence, 
the theory of Gaussian dualities for canonical bosons is 
markedly different and left for a future publication. 

As transformations of statistical mechanics, duali¬ 
ties are valued for mapping strongly-coupled systems to 
weakly coupled ones. Gaussian dualities seem uncon¬ 
ventional from this point of view, since systems of free 
fermions are, by definition, all weakly coupled. To some 
extent, the conceptual mismatch is just a matter of choice 
of language, that is, of physical representation of the sys¬ 
tems under consideration.®® Take for example the sim¬ 
plest model of magnetic ordering, the transverse-field 
Ising chain. The self-duality of this model, a non-local 
transformation of spins, is the prototype of a strong cou¬ 
pling/weak coupling duality transformation. However, 
the Jordan-Wigner mapping transforms the Ising model 
into the Majorana chain, and the self-duality of the Ising 
model into a Gaussian self-duality of the Majorana chain. 

The key property shared by all duality transforma¬ 
tions is symmetry transmutation in face of the locality 
constraint. This phenomenon occurs when a symmetry 
and its dual partner, necessarily a symmetry of the dual 
Hamiltonian, have different physical interpretations. For 
Gaussian dualities in particular we demonstrated trans¬ 
mutation of particle number, fermionic parity, transla¬ 
tion, spin rotation, and time-reversal symmetry in var¬ 
ious space dimensions. Transmutation of particle num¬ 
ber is most conspicuous, since it allows for insulators (or 
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semimetals) and superconductors (with zeroes of the gap 
function) to appear as dual partners. 

Because of symmetry transmutation, Gaussian duali¬ 
ties can establish equivalences of topological insulators 
and superconductors, relating systems classified as in¬ 
equivalent from a single-particle viewpoint. Here we in¬ 
vestigated in detail two paradigmatic examples of such 
pairs of dual partners: the dimerized Peierls/Majorana 
chain (at vanishing chemical potential), and graphene 
which happens to be dual to a weak topological supercon¬ 
ductor. (Transmutation of lattice symmetry is manifest 
in both examples. In particular, the superconductor dual 
to graphene resides on a square lattice.) 

While our list of equivalences is far from exhaustive, 
our approach constitutes a general framework. Hence, 
future research should be focused on searching systemat¬ 
ically for other classes of Gaussian dualities besides the 
ones presented, aiming to obtain all possible equivalences 
across entries in established classification tables of elec¬ 
tronic matter. We can offer two comments as to what 
“possible” should entail. First, we do not necessarilly 
expect equivalences of systems of different space dimen¬ 
sionality, since dimensional reduction by duality is pos¬ 
sible but uncommon.Second, to obtain equivalences 
between topologically trivial and non-trivial systems, it 
will become necessary to obtain the Gaussian analog of a 
holographic symmetry.The Gaussian dualities of this 
paper do not realize holographic symmetries: boundary 
symmetries (zero-energy modes) of a topologically non¬ 
trivial system are mapped to equally localized boundary 
symmetries of its dual partner, and topologically triv¬ 
ial systems are mapped to equally trivial dual partners. 
In short, our dualities preserve the bulk-boundary corre¬ 
spondence (for a useful discussion of this correspondence, 
see for example Ref. [34] and references therein). 

Let us now discuss some possible applications of our 
results. It is interesting to think of dual partner systems, 
with at least one open boundary, as scattering regions 
and ask what happens at the level of dualities if we at¬ 
tach leads to the system. The remarkable answer is that 
it is often possible to extend the Gaussian duality to ap¬ 
ply to the whole system (scattering region plus leads), so 
that its dual partner also has an unambigous interpreta¬ 
tion as a scattering region with leads attached. Depend¬ 
ing on the nature of the leads and Gaussian duality, the 
dual leads may be rearranged in space and/or become 
superconducting. In any case it is now possible to ob¬ 
tain quantitative dual mappings of transport properties 
to further extend the equivalence of insulators and super¬ 
conductors. Notice that it is no problem to add disorder, 
but symmetry transmutation of the disorder ensemble 
should be expected. That is, if the disorder ensemble of 
one system displays some statistical symmetry, the dual 
system will be distributed according to a dual ensemble 
with a dual statistical symmetry of possibly very different 
nature. These brief comments set the ground for inves¬ 
tigating equivalences of statistical topological insulators 
and superconductors.^® 


Topologically non-trivial insulators and superconduc¬ 
tors display zero-energy boundary modes associated to 
the degeneracy of their many-body ground-state energy 
level. In general, it is possible to trace the superconduc¬ 
tor ground degeneracy back to some discrete symmetry 
left over from the breaking of particle conservation. This 
picture does not apply to insulators and so the origin 
of their topological degeneracy is harder to unveil. But 
our equivalence can help here. In terms of its super¬ 
conducting equivalent, the ground degeneracy of a topo¬ 
logical insulator is explained by spontaneous symmetry 
breaking. Moreover, it can be detected and character¬ 
ized numerically by investigating the system in terms of 
Josephson-like physics. 

Finally, we would like to point out a very important 
application of our duality approach to topological mat¬ 
ter. Since our duality transformations are not restricted 
to free fermion/boson systems (as shown in a couple of 
examples of this paper), it is apparent that they may 
become a powerful tool for a potential classification of 
interacting fermion/boson systems. 
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Appendix A: The m-merized Peierls chain 


Let us describe briefly the superconducting model asso¬ 
ciated to the general m-merized Peierls case with m > 2. 
Consider the m = 3 (trimerized) Peierls chain with hop¬ 
pings ti, t 2 , and ts and periodic boundary conditions 
(L = 2M is divisible by m = 3). Then, the dual super¬ 
conducting Hamiltonian is given by 

L 

= -^2 {wj c]cj+i -f Aj c]c]+i -b H.c.) , (Al) 
i=i 


with coupling constants 


ti + ta ti +t2 ^3 + ^2 / . „\ 

Wl = -^- , W2 = -r- , W3 = --- (A2) 


Al 4 = ± 


2 

tl — ts 


A2,5 = ± 


tl — t2 


Aa.e = ± 


ts — t2 


and the rest of the couplings periodically repeat. Repeat¬ 
ing the same line of reasoning, the case m = 4 leads to a 
dual Hamiltonian such as Eq. (Al), where the coupling 
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constants are given by wi = , W 2 = , W 3 = 


, A.-j — 


t3±^2,^^;4 = t3±t4^ = ^^,A3 = 

^ . This clearly shows the difference be¬ 

tween the cases where m is odd from those where m is 
even. In the latter case, the periodicity is always m lat¬ 
tice constants, except for the particular dimerized m = 2 
case where the periodicity is one lattice constant, while 
the periodicity becomes 2 m when m is odd. 


Appendix B: Insulating dual partner of the BCS 
superconductor 

There is an elementary Gaussian duality of an s-wave 
BCS superconductor to a trivial insulator. Let us de¬ 
note by Ck,cr the annihilation fermionic field in momen¬ 
tum space k, and assume Ck = C-k- Then we may define 
the Gaussian duality 

Ck,t ^k.t’ ^-k,i’ (^ 1 ) 

acting non-trivially only on the spin-down fermions. The 
mean-field BCS Hamiltonian is 

H - fj-) c^Ck.a - H ^ (4,t^-k,t + H-C.), (B2) 

k,(7 k 

and gets mapped to 

= '^2 [^kCk,QO-^,/3Ck,/3 ~ (^3) 

k 

with C = I]k(*=k - fJ-), = X, y, z, and 

hx = A, hy = 0, = y,. (B4) 

Particle conservation is restored in at the expense of 
broken (by symmetry transmutation) time-reversal and 
spin-rotation symmetry. 


Appendix C: A Bosonic insulator and its dual 
superfluid 

In recent years, lattice models of bosons have acquired 
new relevance thanks to the spectacular experimental de¬ 
velopment of ultracold atom physics. In this section we 
will focus on a case of practical importance in which only 
three states {\n — I), |h), \n -|- 1)} per lattice site j are 
physically active. Then the particle operators , gj, rij 
act on these states as 

ffj|h-l)=0, 5]|h - 1) = |h), 

g^-|h) = |h - 1), 5]|h) = |h-b I), 

gj\n+l) = y/n\n), g]|n-bl)=0, (Cl) 

and 

nj\n — 1) = (h — l)|h — 1), 
nj\n) = h|h), 

nj|h-b 1) = (h-b l)|h-b 1). (C2) 


These are bosonic operators restricted to a finite(three)- 
dimensional Hilbert space per site. Their algebra and 
application to optical lattices^® has been extensively in¬ 
vestigated in Refs. 37 and 23. Here it will suffice to notice 
that there is a connection to spin S' = 1 operators that, 
in the present paper and because of the duality we will 
apply later on, we take to be 

y|st = S- + iS- = sy n,=S» + fi. (C3) 

Here we will focus on a system of bosonic atoms dis¬ 
tributed with (integer) average particle density h > 1 on 
a one-dimensional optical lattice. We assume that the 
dynamics of the system is described by the Hamiltonian 

H ^ + H.c.) - yUj + V UjUj+i . (C4) 

3 

The chemical potential y fixes h, the interaction H > 0 is 
repulsive, and the total number of particles N = rij 
is conserved. This bosonic t-V model displays a Mott 
insulating phase. 

As a result of Eq. (C3), H has an interesting interpre¬ 
tation as an XXZ spin S' = 1 model in a magnetic field, 

j 

-(/i-2hH)Sj' + ySj'Sj'+i] (C5) 

(up to an additive constant). In the following the chem¬ 
ical potential will be kept fixed at y = 2nV in order to 
eliminate the magnetic field term. 

The quantum phase diagram of the S = 1 XXZ is well 
known.For t = 0, the ground state is antiferromag- 
netically (Ising) ordered in the spin language and in a 
Mott insulating state in the boson language. There is 
a mass gap in the quasiparticle spectrum. According to 
the Haldane gap conjecture, the mass gap remains as the 
Heisenberg antiferromagnetic line t = —2Vjn is reached. 
For the bosonic t-V model this means that the atomic 
system remains in a Mott insulating state. Next we will 
show that the bosonic t-V model is dual to a particle 
non-conserving Hamiltonian with an unconventional su¬ 
perfluid ground state. 

To establish the equivalence of the bosonic t-V model 
to a superffuid, we will exploit a duality transformation 
first investigated in the context of spin S = 1 models, see 
Ref. 39 and references therein. A compact expression for 
the unitary transformation associated to this duality is'*'' 

= (C6) 

j<k 

It follows that 

7y 7yt _ 

(C7) 
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Therefore, the dual spin model =U^HU\ is 


3 


(C 8 ) 


The dual model of bosonic atoms is obtained by rewrit¬ 
ing the spin Hamiltonian in terms of particle opera¬ 
tors. This task is easy once the identity 

-lsy{{s+r + (C9) 

is established. Then it follows that the dual bosonic t-V 


model is described by the Hamiltonian 


H [Y(9]9j+i+g]+i9j) 

-i2^3i9j + 9y)i9j+i + 9jli)nj+i 
The terms 


H^ = 


V 

with fij = rij — n. 


+ 9j ^9jli)nj+i. (Cll) 


(CIO) 


break particle conservation. Bosons are created and an¬ 
nihilated in quartets. Hence the U(l) symmetry of par¬ 
ticle conservation of the t-V model is explicitly broken 
down to a discrete Z 4 symmetry in the dual Hamilto¬ 
nian. The operator = e'^ commutes with , 

and has the interpretation of counting the total number 
of bosons modulo four. 
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